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EXECUTIVE SUMMARY 



This thesis provides closed forms which describe the statistics of the angular 
distribution in azimuth and elevation planes for geometrically based single bounce spheroid 
model. This spatial model can be generated by revolving an ellipse, which lies on the xz 
plane of a 3D Cartesian coordinate system, about the x-axis. This ellipse has semi-major 
axis a and semi-minor axis b , respectively. The base and the mobile stations are located in 
the xy plane with the mobile station being at the center of the spheroid and the base station 
at the origin of the coordinate system. It is assumed that the mobile station is transmitting 
and the base station is receiving. The base and mobile stations are separated by a distance 
D. 

The scattering spheroid model is applicable to a macrocell environment in which the 
received multipath signals at the base station are originated after a scattering from the 
surrounding environment about the mobile station. Assuming that the scatterers are 
uniformly distributed with a constant scatterer density function in the volume V as defined 
by the geometry of the spheroid, we derive the marginal angle of arrival (AOA) probability 
density functions (p.d.f). in the azimuth and elevation plane. 

Using the derived closed forms for both p.d.f.s we generate plots to examine how 
the parameters a, b and D affect the angular distribution of the received waves. From the 
marginal AOA p.d.f. in the azimuth plane we conclude that: 
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1 . The higher the ratio a/D , the lower the probability for the received 
multipath signals at the base station to be confined to a small angular 
region centered about the axis which connects the base and the mobile 
stations. 

2. The angular distribution in the azimuth plane is independent of the length 
of the vertical dimension b of the spheroid . 

3. Comparing the results of the spheroid model for small values of b with 
those of the circular scattering model presented in [Ref. 3], it was 
observed that there exists a significant difference between them. 

Additionally, the effect of the a/D ratio on the r.m.s angle spread in azimuth plane is 
discussed. From the plot results of the r.m.s angle spread in azimuth plane versus a/D, 
we conclude that for values of a/D up to 0.5, there exists a linear relation. However, for 
higher values of a/D, the relationship changes to a nonlinear curve with a positive 
gradient. 



Regarding the marginal AOA p.d.f. in the elevation plane, the following 
conclusions were drawn from the plots: 

1. Similar to the p.d.f in the azimuth plane, the higher the ratio b/D the lower 
the probability for the received multipath components at the base station to 
be restricted in a small angular spread around the axis, which connects the 
base and the mobile stations. 

2. The r.m.s angle spread which also depends on the semi major axis a, has 
a linear dependence with respect to the b/D ratio. 
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I. 



INTRODUCTION 



A. BACKGROUND 

The outdoor radio environment contains a large number of obstacles, such as 
buildings, cars, hills, etc, which obstruct, reflect and scatter the wireless radio signals. As 
a result, the transmitted signal arrives at the mobile station or at the base station from 
many different angles and times of arrival; this is known as multipath. Another common 
phenomenon in cellular systems is the co-channel interference, and it is created because 
in a certain geographic region(cluster) there are other subregions(cells) that utilize the 
same set of frequencies. Multipath fading and co-channel interference can increase 
significantly the bit error rate and they don’t allow high bit rate data services in wireless 
communications systems. The presence of multipath and the motion of the receiver or the 
transmitter causes the received signal envelope to be time-varying and experience a 
Doppler shift that depends upon the angle of arrival of the incoming waves. Concepts 
such as angular spread, time delay spread, and Doppler spread have to be taken in 
consideration during a wireless communication system evaluation. In particular, they 
should be taken into consideration in the design of smart antennas used for diversity, 
beamforming and emitter localization applications. 

The multipath effect alone can be countered by applying antenna spatial diversity 
which implements an antenna array consisting of multi-elements. Spatial diversity is an 
arrangement whereby signals from multiple elements are combined to result in a signal 
that fades less rapidly than the individual element signals. Ideally this can be 
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accomplished if received signals at the two elements are decorrelated. The degree of the 
decorrelation between the two received signals by the two elements can be expressed by 
the spatial cross-correlation coefficient which is a function of the joint p.d.f. 0 (^, 0 ) 

describing the angular distribution of the incoming waves in the three-dimensional space 
at the receiver. 

The interference limited factor can not be overcome by using spatial diversity. 
Switched beam arrays can be used with some success if the angle of arrival of the 
interfering signal is not near the angle of arrival of the desired signal [Ref. 1], Only the 
beamforming arrays have the potential to completely suppress the co-channel interference 
( if the number of elements exceeds the number of co-channel interferes) improving the 
wireless communication system performance. The performance of these type of smart 
antennas depends on the spatial characteristics of the wireless radio channel. The term 
“spatial characteristics of the radio channel” is used to mean the statistics of the received 
signal envelope in the spatial domain. These spatial characteristics may not always be 
available from measurements. In such a case one needs to devise geometric models to 
depict the scattering process in the spatial domain. Another area where knowledge of 
angular distribution is valuable is in the determination of the power spectral density [Ref. 
1]. Therefore the angular distribution of the received multipath components is important 
in evaluating the performance of radiowave link between the transmitter and the receiver 
in wireless communications. 
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Spatial channel models have been created and applied in order to be able to derive 
the statistics of the received signal envelope for both macrocell and microcell 
environments. The term macrocell enviroment assumes that the scatterers surrounding 
the mobile station have the same height or are higher than the antenna of the mobile. 
This results in a received signal at the mobile station which arrives from all directions 
after bouncing from the surrounding scatterers. In contrast, the base station’s antenna is 
located higher than the surrounding scatterers. Therefore the received multipath signals 
at the base station originate after scattering from the area near to the mobile station 
[Ref.2]. In microcell enviroment both base and mobile station are deployed in the same 
height as the surrounding scatterers. This implies that the scattering process takes place 
near to the base and mobile station. The characteristic difference between the types of 
environments is that the angular spread of the received signals at the base station is 
confined to a small region for macrocell environments. In the microcell enviroment this 
region tends to be larger. 

In Table 1, which is a modified version of a table from [Ref 2.], are listed some 
representative spatial channel models. 
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SPATIAL CHANNEL 

MODEL 


DESCRIPTION OF THE 
MODEL 


APPLICATIONS 


1. Lee’s Model 


Scatterers are evenly spaced 
on a circular ring about the 
mobile at a distance R. Each 
of the scatterrers represents 
the effect of many other 
scatterers and hence is 
referred to as an effective 
scatterer. The base station is 
deployed a distance D from 
the mobile station. 


Predicts correlation 

coefficient using a discrete 
angle of arrival(AOA) 
model. A modified version 
of this model can account 
for a Doppler shift, 
assuming that the scatterers 
are moving with an angular 
velocity. 


2. Discrete Uniform 

Distribution 


A model similar to Lee’s 
model with the N scatterers 
evenly spaced over an 
AOA range. 


Predicts correlation 

coefficient using a discrete 
AOA model. The 

correlation predicted by this 
model falls off more quickly 
than the correlation in Lee’s 
model. 


3. Geometrically Based 
Circular Model 


In this model it is assumed 
that the scatterers are 
uniformly deployed within a 
radius R about the mobile. 
The locations of the 
scatterers are defined by a 
spatial scatterer density 
function. 


It is suitable for macrocell 
environments and predicts 
the joint time of 

arrival(TOA) and AOA and 
marginal density functions 
at both base station and 
mobile station. 


4. Geommetrically Based 
Elliptical Model 


Scatterers are uniformly 
distributed in an ellipse 
where the base station and 
mobile station are the focii 
of the ellipse. 


It can be determined the 
marginal AOA, TOA, the 
joint TOA and AOA density 
function. It is applied in 
microcell enviroment. 


5. Gaussian Wide Sense 
Stationery Uncorrellated 

Scattering (GWSSUS) 


N scatterers are grouped 
into clusters in space such 
that the delay differences 
within each cluster are not 
resolvable within the 

transmission signal BW. 


Provides an analytical 

model for the array 

covariance matrix. 
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6. Gaussian Angle of 
Arrival (GAA) 



Special case of the 
GWSSUS model with a 
single cluster and angle of 
arrival statistics assumed to 
be Gaussian distributed 
about some nominal angle. 



Provides an analytical 
model for the array 
covariance matrix. 



Table 1. Description of some representative spatial channel models. “From Ref. [2].” 

Figure 1 through 4 demonstrate the first four spatial channel models, listed in 
Table 1. 
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Figure 2. Discrete Uniform Distribution Model “From Ref. [2]”. 
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Every spatial channel listed in Table 1 has one common feature; the received 
signal at the base station is assumed to arrive from a two dimensional plane connecting 
the tips of the transmitter and the receiver. Hence the angle of arrival (AOA) includes 
only the azimuth angle information. 

In this thesis we consider a 3D spheroid model as shown in Figure 5 This 
spheroid can be generated by revolving an ellipse, which lies on xz plane, with semi- 
major axis and semi-minor axis a and b, respectively, about the x-axis. The base station 
and the mobile station are located in the xy plane with the mobile station at the center of 
the spheroid, and the base station at the origin. It is assumed that the mobile station is 
transmitting and the base station is receiving. The distance between the two terminals is 
D with the requirement that the semi-major axis of the spheroid be less than or equal to 
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the distance D. In the model considered here it is obvious that the angle of arrival (AOA) 
of the received signals includes both the azimuth and elevation plane coordinates. 




Figure 5. The spheroid geometry. 
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B. OBJECTIVE 

Assuming that the scatterers are uniformly distributed around the mobile station 
as defined by the geometry of the spheroid, the objective is to derive the marginal 
probability density function (p.d.f.) angles of arrival (AOA) in the azimuth and in the 
elevation planes at the base station. After extracting these quantities, we investigate how 
the angular spread in azimuth and elevation plane are related with the dimensions of the 
spheroid and the distance between the transmitter and the receiver, as well as how they 
are correlated with each other. 
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II. PROBLEM FORMULATION 



A. INTRODUCTION 

In this chapter we consider a single bounce scattering spheroid model. The model 
can be used in a macrocell environment where the base station antenna height is 
relatively high. As a result there will be no scattering from scatterers near the base 
station. Figure 6, shows the spheroid scatterer density geometry and the notation used to 
derive the marginal p.d.f.s. As is shown in the Figure 6, the location of each scatterer is 
defined by the Cartesian coordinates (x,y,z) or the spherical coordinates (r b ,<p b ,9 b ), or 

(r b ,(f> b ,j3), where (3 represents the complementary angle of 9 b in elevation plane. 
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Figure 6. The spheroid scatterer density geometry. 

The method for the derivation of the joint TOA/AOA p.d.f. and marginal AOA p.d.f. for 
both azimuth and elevation planes is based on the approach presented in [Ref. 3]. In the 
Figure 6, it is assumed that scatterers are distributed randomly in the volume V of the 
spheroid according to the spatial scatterer density function f x y z (x,y,z). It is further 

assumed that the waves arrive directly at the base station after undergoing scattering only 
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once from within the volume V. Multiple scattering of waves is not considered in this 
thesis. 

Before proceeding with the determination of the p.d.f.s in the azimuth and the 
elevation plane we would like to present the formulation for the joint TOA/AOA p.d.f 
and marginal AOA p.d.f in the general case. In section B we start with the derivation of 
the joint TOA/AOA p.d.f.. In section C, we show formulation for the marginal AOA 
p.d.f.s. In the following sections D through G we derive the marginal AOA p.d.f in 
azimuth and elevation plane for the spheroidal model. 

B. DERIVATION OF THE JOINT TOA/AOA P.D.F.’S 

Figure 7, shows a plane wave impinging on a base station after scattering from an 

obstacle, represented by a point ‘S’ in space with the ordered rectangular triples (x,y,z). 
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The scatterer probability density function f x y z (x,y,z) can be also expressed 

terms of spherical coordinates systems {r b ,(f) b ,6 b ) as follows. The equations relating the 
spherical coordinates to the Cartesian coordinates are: 

( 1 ) 



r b - a /* 2 + y 2 + z 2 



6 b - cos 



-1 



■Jx~+y~+z~ 



( 2 ) 
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= tan 1 (— ) 

v jr' 



( 3 ) 



x-r b sin 6 b cos (j) b 
y - r b sin 9 b sin 



Z = r b cos 0 b 



(4) 



(5) 



( 6 ) 



where (x, y, z) and ( r b ,(j) b ,d b ) denote the location of the scatterer in Cartesian and 
spherical coordinates respectively. The joint p.d.f f ^ e {r b ,<f> b ,9 b ) can be found using 
the well known transformation [Ref. 4]: 



x = r b sin 9 b cos (j> b 

fr h ,<l> b ,oS r b'fo' e b)= J ( r b’ < t>b’ d b)fx,y,z( X 'y' Z )y = r b sin 9 b Sin ^ 

z = r b cos 9 b 

where j(r b ,</> b ,9 b ) is the determinant of the Jacobian transformation given by: 



(7) 



j{ r b’0b'&b)~ 



dx dx dx 



dr b d<p b dO b 
dv dv dv 



dr b d<p b d9 b 

<t <k <k 



dr b d<p b dd b 



= r b sin 9 b , 0 <9 b <n 



( 8 ) 



Substituting the equation ( 8 ) into (7) we get for the scatterer probability density function 
/ 0 e (r b ,<P b ,0 b ) with respect to spherical coordinates: 



AaA ( r b ,<PbA) = r b sin e bfx,y,z ( r b sin 6 b cos <t>b . r b sin e b sin h > r b cos e b ) • 



(9) 



If we want to express the scatterer probability density function f r ^ e (t~ b ,(j) b ,0 b ) 



with respect to (Z,<p b ,6 b ), where T denotes the time delay of the multipath component, 
we have first to find the relationship between r b and the T. Applying the law of cosines 
to the triangle BS,S,MS we get the following equation: 
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r? = r b + D 2 - 2r b DcosS 



( 10 ) 



where the cos 8 is called the ‘cosine direction’ given by: 

cos 8 = sin 6 b cos <p b . 

The delay corresponding to the total path traveled by the wave is: 
= 7 +V r i 2 + -D 2 + 2r b Ds\nd b cos(f ) b ) 



( 11 ) 






(12) 



where ‘c’ is the speed of propagation of waves. 

Squaring both sides of the equation (12) and solving for r b gives: 



D 2 -t 2 c 2 



(13) 



Vb 2(3 sin W b cos ~ij> b -rc) 

The joint TOA/AOA p.d.f. f T ^ 0 ("T, <p b ,0 b ) can be found using the Jacobian 



transformation [Ref. 4] according to: 

f (-a a fn,.f b .e b ( r bA^b') 

} ’M > - \H > iAAH 



2 ° 

P-TC 



r b 2(Dsin 6 b cos <p b -rc) 



(14) 



where \j(r b ,0 b ,6 b ] is the determinant of the Jacobian tranformation given by: 



| j{ r b^b^b \ - 



2a. 

dr 



-1 



2(Dsing 6 cos 0 b -Tc) 2 
D 2 c+z ,2 c 3 -2 '.tc 2 D sin & b cos $ 



(15) 



Substituting the equation (15) into (14) we get for the joint TOA/AOA p.d.f. 

fz A .e, 0> A) 

f ( T A o \ _ D 2 c+t 2 c 2 ~2rc 2 Ds,me b cos<p b f [ D 2 -t 2 c 2 a (u 

Jr,<t> b ,e b \ T ’0b’ b) 2(Dsin& b cos<t> b -Tc) 2 Tb ’^ b '° b ' 2 ( Dsin &b cos J b -re) ’ rb ’ *>)■ 

Using the last equation (16) and the equations (9) and (13) we can express the joint 

TOA/AOA p.d.f. / T p e * n terms of original scatterer probability density 



function f x _(a, ;y,z) as follows: 
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r , , n \ n (D 2 -T 2 c 2 ) 2 (D 2 c+T 2 c 3 -2rc 2 DsinO b cos 0 b ) r / „ , _ _ \ 

fr.9 b .e h ^4bA) = sm e b — — 71 fx,y,z\ r b sm &b cos (/> b , r* sin d b sin <p b , r b cos 6 b ) • 

° ° b(Z 9 sin & b cos <p b -zc) 



(17) 

Note that the equation (17) expresses the joint TOA/AOA p.d.f. f r ^ e ( T,<p b ,d b ) as 
observed at the base station in terms of an arbitary scatterer probability density function 

fx,y.z( X 'y* Z )- 

When the scatterers are uniformly distributed within an arbitrarily space A with a 
volume V, then the scatterer probability density function f xyz (x, y, z) is given by: 



, / \ , },x,y,zeA 

fx, y ,z\ x >y ^)= il , 

0, otherwise 



(18) 



If we consider this special case, the joint TOA/AOA p.d.f. f 9 (r,<p b ,6 b ) can be 



written: 



fr,<p b M b ( T ^b^b) = sin0 b 



(/> 2 -r 2 c 2 ) 2 (P 2 c+r 2 c 3 -2t<: 2 Z)sin^,cos^,) 
8V (d sin 6 b cos <p b -zc) 4 



where the range for the values of T,0 b ,O b has to be defined. 



(19) 



C. DERIVATION OF MARGINAL AOA PROBABILITY DENSITY 
FUNCTION 

The marginal AOA p.d.f. with respect to spherical angle coordinates could be 
found by integrating the joint p.d.f. / ^ ejjb^b^b) » as gi ven in equation (9) with 

respect to r b over the range {(p b , d b ) to r b {<p b , 9 b ) , and then either with respect to 
eM over the range 6 b (<p b ) t0 °b 2 (<Pb ) or with respect 0 b (6 b ) over the range 
\{ e b) to 0 b ^ {6 b ) depending on whether the marginal A.O.A p.d.f. is to be found in 
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azimuth or elevation plane respectively. The terms r b {(f) b ,G b ) and r b (<p b ,G b ) are the 
lower and the upper limits of integration and depend on the azimuth angle ^ and the 



elevation angle 6 b . The values of these limits can be found from the intersections of a 
straight line 6 b = constant, (j> b = constant with the scatterer volume. The limits of 
integration in 6 b , G bx and.G bl for the marginal AOA p.d.f. in azimuth plane can be 
defined as the smallest and the largest 6 b -values that bound a region R, which is the 
intersection between a plane at a specific (f> b value and the scatterer region volume. The 
above statement can be expressed in mathematical expressions as follows. 

The marginal AOA p.d.f. {ff> b ) in azimuth plane is given by: 



0b 2 i^b) r b2^b^b) 

J J rg sin e b f Xt y tZ (r b sinG b cos (p b ,r b sin G b sin (j) b , r b cos G b )dr b d<p b . (21) 

<%1 (&b ) % (0b A > ) 

The limits for (j) b integration, which depend on the elevation angle G b , can be 
determined from the smallest and the largest azimuth angles that bound the region of 
intersection between the cone G b = constant and the scatterer volume. 




^2 (0b) r lri(0b^b) 




The marginal AOA p.d.f. f 0 (d b ) in elevation plane is given by: 
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Equations (20) and (21) express the AOA marginal p.d.f.’s for any scatterer 



probability density function f x >z (*,y,z). When the scatterers are uniformly 
distributed in a volume with a constant probability density function f x _(jc, y,z), the 
marginal AOA p.d.f. (^) is given by: 



fyn (&b) r b-> (fib fib) 



f<p h (fib ) “ V 



'fib 



J r~ sin 0 b dr b d0 b = 

Pb) r b\ ifafib) 




0 h* e b ) “ r l (<t>b’ e b )l sin d b de b 



( 22 ) 

and the marginal AOA p.d.f. f Q ( O b ) is given by: 

fib 2 (&b ) r b 2 (fib &b ) fib^ (&J > ) 

fd b (°b) = V J J r* sin 6 b dr b d(j) b = ^ J [/£ A ) “ (A , 6 b )Jsin 6 b d<p b . 

<Pb l ( e b) r b X (fib fib) <Pbi i°b ) 



(23) 

D. MARGINAL AOA P.D.F. IN AZIMUTH PLANE FOR THE SPHEROID 
SPATIAL CHANNEL MODEL 

Before we proceed to find the marginal AOA p.d.f. in the azimuth plane, we have 
to define the equation which describes the geometry of the spheroid in Cartesian and 
spherical coordinates. The spheroid in Cartesian coordinates is given by: 



j^Dhl | .- 2 _ i 

a 2 b 2 ~ 

and in spherical coordinates is given by: 

(a 2 cos 2 6 b + b 2 sin 2 6 b )r 2 - 2 Db 2 r b sin d b cos</) b + b 2 [p 2 -a 2 ) = 0 . 



(24) 



(25) 
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